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We study the filling of states in a pure hopping boson model on the comb lattice, a low dimensional 
discrete structure where geometrical inhomogeneity induces Bose-Einstein condensation (BEC) at 
finite temperature. By a careful analysis of the thermodynamic limit on combs we show that, unlike 
the standard lattice case, BEC is characterized by a macroscopic occupation of a finite number of 
^SJ ' states with energy belonging to a small neighborhood of the ground state energy. Such remarkable 

\ feature gives rise to an anomalous behaviour in the large distance two-point correlation functions. 

■ Finally, we prove a general theorem providing the conditions for the pure hopping model to exhibit 

the standard behaviour, i.e. to present a macroscopic occupation of the ground state only. 

D ■ 

C/^ ■ I. INTRODUCTION 
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The most recent experiments on Bose-Einstein condensation (BEC) stimulated a large wealth of theoretical work 
aimed at a better understanding of the basic properties of such interesting phenomenon [1]. In particular the possibility 
of confining ultra-cold bosonic atoms within optical lattices, along with the most recent studies on arrays of Joscphson 
junctions, arouse a great interest in bosonic models defined on Euchdean lattices [2,3]. 

■ Even more interesting from this point of view would be the arrangement of the Josephson junctions, or possibly 
, optical traps, into more complex networks [4]. Indeed it was recently put into evidence [5-7] that, due to topological 

inhomogeneity, BEC at finite temperature can occur on low dimensional structures, such as the comb lattice [8], even 
04 , in the absence of external potential. Such results were obtained mainly in the thermodynamic limit (i.e. for structures 

■ of infinite size), but it is clear that a deeper understanding of the onset of the phenomenon on finite-size structures 
is needed. This is especially true in view of the experimental research which is currently being developed in this field 
[4]. 

One of the distinctive features of usual BEC (in continuous Euclidean geometry, possibly in the presence of harmonic 
potentials) is the macroscopic occupation of a single quantum state. More precisely it is possible to prove that the 
filling of any excited state vanishes in the thermodynamic limit. As for general networks, Refs. [5-7] mainly deal 
Q with the thermodynamic aspect of the problem. In particular it is shown that on a low dimensional inhomogeneous 
O ■ network, such as the comb lattice, BEC is characterized by the macroscopic occupation of the states belonging to 
\ an arbitrarily small energy neighborhood of the ground state. This is a more general condition, since it does not 
necessarily entail that only the ground state features a macroscopic occupation. 

In the following we indeed show that the macroscopic occupation of the ground state does not entirely account for 
^ . BEC on the comb lattice. Even in the thermodynamic limit, the filling of the structure is completely described only 
■ - ' if a finite number of states belonging to a small neighborhood of the ground state is considered. This result is proven 
by exactly evaluating the occupation of the quantum states of mesoscopic comb lattices, with large but finite size. 
The spectrum of such structures is nearly continuous, but the energy levels are still distinguishable since they pertain 
to orthogonal wave functions which differ on a macroscopic scale. 

The paper is organized as follows. In Section II we introduce the pure hopping model [5] which describes the low 
coupling limit for the Josephson junction arrays or a model of non interacting atoms confined in an optical lattice. 
Using the Van Hove spheres, we define the thermodynamic limit for a general infinite discrete structure. Then we 
recall the definition of the low energy hidden spectrum and we give the general conditions on the graph spectra for the 
condensation in an arbitrarily small energy region [7] . In Section III we consider the problem of the states filling for 
the pure hopping bosonic model comparing the known results for the 3-dimensional lattice with the numerical results 
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obtained on the comb graph. We show that in the first case only the lowest energy state presents a macroscopic 
occiipation, while the comb graph exhibits a different behaviour since a macroscopic occupation of many states is 
present. In Section IV we solve analytically the problem of the states filling on the comb graph proving that, in the 
thermodynamic limit, also the first excited states are filled by a finite fraction of particles. Furthermore we show that 
the numerical data perfectly fit the analytical predictions. In Section V we prove that the macroscopic filling of many 
states gives rise to an anomalous behaviour of the large distance two point correlation functions. In Section VI we 
give the spectral conditions for the pure hopping model to present a macroscopical occupation of a single quantum 
state and in Section VII we present our conclusions. 



II. THE BOSONIC HOPPING MODELS ON GRAPHS 



Graph theory provides the most natural mathematical description of generic discrete networks [9]. A graph is a 
countable set V of sites i connected pairwise by a set E of unoriented links = Two sites joined by a 

link are called nearest neighbours. The topology of a graph is fully described by its adjacency matrix Aij, with 
Aij = 1 if (i, j) is a link of the graph and Aij = otherwise. A walk in G is a sequence of concatenated links 
{{i, k), {k,h) . . . , (n, m), (m, k)} and a graph is said to be connected if for any two sites there is always a walk joining 
them. The length of a walk is the number of links appearing in the sequence, and the length of the shortest walk 
joining two sites is called the chemical distance between them. The latter defines the intrinsic metric every connected 
graph is endowed with. The Van Hove sphere Sr.o of center o and radius r is the set of sites whose chemical distance 
from site o is equal or less than r. In the following we will consider only connected graphs with polynomial growth, 
where the number Nr^o of sites within Sr^ grows at most as a power of the radius. This choice ensures that the 
discrete structure can be embedded in a finite dimensional Euclidean space. Since for this class of graphs it is possible 
to prove that the thermodynamic limit is independent from the choice of the center of the sphere [10], in the following 
we will drop the relevant subscript. 

Let us recall some basic results obtained for the pure hopping model on a generic graph (a complete description of 
the subject can be found in [6,7]). On a graph G the pure hopping model for bosons is defined by the Hamiltonian 
[5]: 

where and are the creation and annihilation operator at site i ([a^, aj] = Sij) and Aij is the adjacency matrix of 
the graph. 

To study the behaviour of (1) in the thermodynamic limit, where BEG can arise, we restrict the model to and 
analyze its properties as the radius of the Van Hove Sphere goes to infinity. The model restricted to Sj- is defined by 
the Hamiltonian 

W = Y,-tA\^a\aj (2) 

where A^^ = Aij if i,j G Sr and A^j = otherwise. The corresponding normalized density of states p^{E) is 

P'{E) = ^Y.^{E-El) (3) 
k 

where E^. are the eigenvalues of —tA\j and is the number of sites within the sphere. The function p{E) is defined 
to be the thermodynamic density of states of —tAij if it satisfies the following condition: 

lim [ \p' (E) - p{E)\dE = 0. (4) 

r— *oo J 

In general the asymptotic behavior at low energies of the thermodynamic density of states is described by a power 
law of the form 

p{E) {E - Em)"^-^ for E^ Em (5) 
where E^, = Inf(Supp(/9(S))) and a is an exponent determined by the topology of the graph [11,12]. 
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On inhomogcncous structures the density of states can present some interesting anomalies, such as the hidden region 
of the spectrum defined in [6] and [7] . A hidden region of the spectrum consists of an energy interval [Ei , E2] such 
that [Ei,E2] n Supp{p{E)) = and limr^oo -^[Si e^] ^ ^' where N^^i e^] number of eigenvalues of —tA\- in 

the interval [Ei , E2] ■ Notice that in general N^^^ can diverge for r ^ 00 and the eigenvalues can become dense 
in [El, E2] in the thermodynamic limit. Therefore the presence of a hidden spectrum is a far more general property 
than the existence of a discrete spectrum, where a finite number of states fills the spectral region. An interesting 
example of this new kind of behaviour is exhibited by the the comb lattice [5,6] (sec figure 1) which will be studied 
in detail in the following sections. We now define the lowest energy level for the sequence of densities Pr{E), setting 
Eq = Inffe(£^) and Eq = limr->oo Eq- In general, £^0 < Em- If Eq < Em, then [Eo,Em] is a hidden region of the 
spectrum which will be called low energy hidden spectrum. 



FIG. 1. The comb graph, which is obtained connecting to each site of a linear chain, called backbone, a 1-dimensional chain 
called finger. 



The correct way of taking the thermodynamic limit consists in adjusting the population of the system so that 
the filling N/Nj. is set to a fixed value / as the radius of the Van Hove sphere goes to infinity. In particular, in the 
macro-canonical ensemble, the equation that determines the fugacity 2; as a function of /? = T~^, f and r is: 

/ = J ,-ijE_^ P'^(E)dE. (6) 

Setting Eq = yields < z{f, f3,r) < 1. A system presents BEC at finite temperature if there exists a temperature 
Tc such that for all T<Tc lim^^oo z{f, /3, r) = 1. 

The general conditions for the occurrence of BEC at finite temperature is strictly related to the properties of p{E) 
[7] . In particular it is proven that BEC arises cither in models presenting a low energy hidden spectrum or in models 
where the parameter a appearing in equation (5) is larger than two. The comb lattice is an example of the former 
situation, whereas the latter condition is satisfied by any Euclidean lattice with dimension larger than two, since for 
such a structure the parameter a coincides with the Euclidean dimension. 



III. THE SINGLE STATE FILLING 



BEC on Euclidean lattices is characterized by the non-analytic behaviour of the thermodynamic functions at the 
critical temperature and by a macroscopic filling of a single quantum state. The definition of BEC we gave in the 
previous section is strictly connected with the singular properties at Tc- For example in the thermodynamic limit 
the average energy per particle is: 

= 7 / ,-ijE_^ PiE)dE (7) 

where z = 1 for temperature below the critical temperature, whereas for T > Tc it is a function of T and / determined 
by (6). Therefore the equation (7) presents the typical non-analytic behaviour of the thermodynamic functions of 
BEC at the critical point. 
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Let lis now consider the filling of the states. Below Tq any arbitrarily small energy region around is filled by a 
finite fraction of particles, as it is known by classical results on BEC. More precisely, for all e > the fraction of 
particles with energy in the region [Eq, Eq + e] is always greater than the positive quantity Utd- 

^TD = f - j ^^p{E)dE (8) 

where £^0=0 and the subscript TD stands for thermodynamic. 

Since Eq is the only real number belonging to [Eq, Eq + e] for any value of e, one would argue that the only possible 
way to satisfy the condition rif > jitd > Ve > is to fill the lowest energy state with a finite fraction of particles. 
However if we look carefully to how the thermodynamic limit is performed, the question is not so trivial. Indeed, 
if the model does not present a gap (this is also the case of the usual condensation on lattices), the energy of the 
first excited states tends to Eq in the thermodynamic limit and it belongs to [Eo,Eo + e] for all values of e. For 
instance, in order to determine the filling ni of the first excited state Ei, one first has to solve equation (6) for z{r) 
and then evaluate rii = limr^oo exp(/3i^i(r)) — The result of this limit depends on how oo, 

Ei{r) Eq = and ^(r) — > 1 in the thermodynamic limit. 

This limit is trivial on lattices: there the pure hopping model does not present a gap and it is known that below 
Tc the lowest energy eigenstate is the only one with a macroscopic occupation [13]. Let us first consider the case of 
the three-dimensional lattice. The spectrum of a finite lattice of = Nj- sites is given by: 

ai = <6t-2t cos — 2t cos — 2t cos — ^ \ (9) 

the equation determining the fugacity z for each finite lattice is: 

.1-2 1 ^ — r Z , , 

^^IfiY^z^l? ^ e^^^ - z ^ ^ 

where no = L~^{z~^ — is the filling of the ground state and L~^(z~^ cxp{/3Ek) — is the filling of a state of 
energy Ei~. Below the critical temperature equation (10) has a solution in the thermodynamic limit only ii z 1 for 
L — > oo. In particular one has that: 

z~l— forL^oo (11) 
Substituting (11) in (10) one obtains no = and 



"0 = /-/ ^^^PiE)dE- 



(12) 



Hence no = nxo ^mcl the lowest energy state is the only macroscopically filled state. As a further check of this result 
we can explicitly evaluate the filling ni of the first excited state in the thermodynamic limit. According to equation 
(9) the energy of the first excited state is Ei = t{2 — 2cos(27r/L)) and, from equation (11), its filling is 

1 z 11 
m = lim —- — -= lim -7? -,--„ = (13) 

in the thermodynamic limit. 

In figure 2 we plotted utd (obtained from (8) with / = 1), no (the filling of the lowest energy state), and n^ for 
two different values of e. To determine no and n^ we first evaluated the exact spectrum of a finite lattice consisting of 
Nr = 125000 sites, then we obtained z{l,/3,r) by numerical inversion of equation (6). This allowed us to evaluate the 
filling Hfc = N~^{z~^e^^'' — 1)^^ of the energy level Ek- From figure 2 we have that the differences between the plots 
are very small and they can be ascribed to finite size effects. These numerical results confirm the known property of 
BEC on regular lattices of presenting a macroscopic occupation only in the state of lowest energy. 
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FIG. 2. The state filling in the 3-dimensional lattice as a function of the temperature measured in units of t. The dotted line 
represents the theoretical value of utd obtained from (8). The solid line is a numerical evaluation of the filling of the ground 
state no. The dashed and dashed-dotted lines are a numerical evaluation rie for two different small values of e: e = 0.02i, 
e = 0.05t. Numerical data are obtained from a lattices of Nr = 125000 sites. 



Let us focus on the case of the comb graph proving that on inhomogeneous structures, due to the presence of hidden 
regions in the spectrum, a macroscopic fiUing of states of energy arbitrarily near to the ground state is possible. In 
figure 3 TiTD, and for a comb graph consisting of 40000 sites are shown. In this case it is evident that, since the 
fiUing of the ground state (the curve no) is lower than titd, there must be macroscopically filled states other than 
Eq. Prom the general result (8) the energies of these states must tend to Eq as L ^ oo. Actually, the curves 
representing the filling relevant to the energy region [£^0, -£'0 + e] is very close to titd even for small values of e. 




FIG. 3. The state filling in the comb lattice as a function of the temperature measured in units of t. The dotted line 
represents the theoretical value of htd obtained from (8). The solid line is a numerical evaluation of the filling of the ground 
state no. The dashed and daslicd-dottod linos arc a numerical evaluation Ue for two different small values of e: e = 0.02t, 
e = 0.05t. Numerical data are obtained from a comb of A',. = 40000 sites. 



IV. STATES FILLING ON THE COMB GRAPH 



On the comb graph, the anomalous behaviour of the filling can be explained by carefully considering equation (6), 
which determines the fugacity for each finite-size comb. The energy spectrum of the pure hopping model on a finite 
comb graph consisting of Nj. = L x L is the union of four sets, ctp = {Eq} U a- U ao U cr+ [6]. Setting Eq = 0: 

a- = lEk-- 
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ao = {Ek=t(V&-2cos^)] (14) 

I V ^ J ) k=l,2,L-l 



0-+ 

The degeneration d{Ek) of the states of energy Ek G ct- and Ek € cr+ is 2, while for the states of energy Ek G ctq 
we have d{Ek) = L. Hence, in the thermodynamic Umit, C7_ and cr+ are fiUed by a vanishing fraction of states 
{2{L — 1)/(4L^)) and they belong to the hidden spectrum [6], whereas gives rise to the spectral region of measure 
one, since the relevant fraction of states is L{L — l)/LF'. ctq reproduce the spectral density of a linear chain and 
its states are completetely delocalized. On the other hand the states in cr_ and a+, which characterize the typical 
behaviour of the comb graph, are localized along the backbone presenting an exponential decay in the direction of 
the finger. 

The equation determining the fugacity z for a finite comb is 

/ = riQiz, L) + n„_{z, L) + n„„{z, L) + n„^{z, L) (15) 

with < z < 1. Here no{z,L), na^{z,L) , na_{z,L) and nao{z,L) represent respectively the fraction of particles in 
the ground state and in the three spectral regions. More precisely 

Mz,L) = ^^d{E)-^^ (16) 

Eea- 

where d{E) denotes the degeneracy of the energy state E. Let us analyze these four contributions separately. 
In the thermodynamic limit the filling of the ground state, 

1 z 

noiz) = lim no{z,L)= lim — , (17) 

L— >(X> i— >(X> h 1 — Z 

and the filling of the low hidden region, 

n,_(z)= lim n<,_(z,L)= lim 4 V U" V*(^-2 V^W^) _ 1 ^ (Ig) 

n=l 

are different from zero if and only if limj^^oo z = \. The filling of the spectral region of measure 1 is given by: 

1 ^ z 1 'Y /" 1 

n., (.) = ^lim n., (., L) = £m^ - g ) _ ^ = ^ ^^ ^-1^/^^78-2 cos fe)) _ i = / .-^e^B _^ pWE 

(19) 

na^iz) is a finite positive number for each value of z, < 2; < 1. For the high hidden region o-+ we obtain: 

2 z 2 z 
Q<naAz)= lim naAz,L)= lim — , < lim — — ^ =0 (20) 

- > L-,00 + ^ ' ' L-,00 L2 g/3t(V^+2 ^l+cos2 2p) _ ^ L-,00 L e/3«(v/8+2) -z 

Hence {z, L) can be neglected. 

Let us consider the behaviour of equation (15) in the thermodynamic limit. Above the critical temperature it is 
satisfied by a value of z = z' < 1 so that, according to equations (17) and (18), no{z') = ncr_ {z') = 0. On the other 
hand for T < Tc equation (15) can be solved only letting 2; — > 1 when i — > 00. In particular we obtain a solution of 
equation (15) if and only if: 

z~l--^ forL^oo (21) 
Let us study in details the behaviour of no{z) and n<j_ (z) when z tends to 1 as in equation (21). We get 

no = ^ (22) 
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and the fraction of particles in cr_ can be exactly evaluated by summing the corresponding series: 



hm -2 



(Z,-l)/4 

E 



c 



i2 



y ^ 



csch a [a cosh a — sinh a] 



where 




Then below Tc equation (15) for z becomes an equation for the new variable 5, 

+ cscha [acosha — sinha] = / — 71^0(1) = utd 



(23) 



(24) 



(25) 



which can be numerically solved. Hence, in order to obtain equation (25) we had to evaluate the thermodynamic 
limit of equations (17), (18), (19), (20). We remark that only in the case of the spectral region of measure one it is 
possible to replace the sum by an integral, whereas in the other cases the sums have to be explicitly calculated. 



Equation (25) explains the main properties of BEC on the comb lattice. First of all in this case utd = f 



,(1) 



is the sum of two contributions, the first due to the particles in the lowest energy state uq = 6~ , the other due to 
the particles in the low hidden region, n^-. = cscha [acosha — sinha]. Hence, on the comb graph, the filling of 
the ground state is smaller than utd, as it has been numerically shown in figure 3. 

A second point is that only states with arbitrary small energy contribute to htd- Indeed the result of equation 
(23) does not change if we force the index n to be lower than eL, where c is au jubitrarily small fixed parameter. 
This means that, Ve > 0, only states with energy smaller than t(\/8 — 2 -^/l + cos2(27re)) contribute to utd- Finally 
substituting S^^ with no in equation (25), we obtain an exact relation between no and utd- In figure 4 we checked 
that the numerical results for the comb graph presented in the previous section are perfectly reproduced by the exact 
calculation. The dashed-dotted line is obtained by first evaluating ng and then using equation (25) to get nTnii^o)- 
These data are very close both to (filling of a small energy region around Eq) and to the theoretical value of utd- 




FIG. 4. Check of equation (25). The dotted line represents the theoretical value of utd obtained from (8). The solid line is 
a numerical evaluation of the filling of the ground state no- The dashed line is a numerical evaluation for e = 0.02t. The 
dash-dotted line is obtained evaluating utd from (25), where 6 and a are given by (24), (22) and no is given by the numerical 
data of the solid line. The temperature is measured in units of t and numerical data are obtained from a comb of Nr = 40000 
sites. 



V. COHERENCE PROPERTIES OF THE CONDENSATE 



The classical Bosc-Einstcin condensates on homogeneous lattices present important large scale coherence properties 
due to the macroscopic filling of a single quantum state. One of the most effective ways to put into evidence this 
relevant feature of BEC is the study of the long distance correlation function C{i,j) defined by: 
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k 



1 



(26) 



Where tpkU) is the eigenvector of —tA^^ (2) of eigenvalue E^.. In particular, considering the correlation function 
between two sites at a macroscopic distance (i.e. rij ~ r, r is the radius of the Van Hove sphere), below the critical 
temperature in the thermodynamic limit we have: 



C{i,j) ^ N^nTDi^l''{i)Vo{j) for r ^ oo 



(27) 



Equation (27) shows that below Tc there arc NrUTD particles in the ground state On regular lattices, the 

wave function of the ground state is constant and the correlation function at large distances does not depend on rij . 
In particular for the three dimensional case {Nr = L^) C{i,j) ~ {L^nTD){L^)^^', {L^nTo) is the asymptotic number 
of particles of the condensate and {L^)~^ is the normalizing factor given by the wave function V'o(i)- 

The large scale coherence properties of the Bose-Einstein condensate on a comb graph can be studied by calculating 
the correlation functions between the sites i and j of the backbone at a chemical distance — dL, (this way if 
the size of the macroscopic system is defined to be one, we are considering sites at distance d). Let us evaluate the 
contribution to (26) of each single spectral region by using the exact wave functions ipj^{j) of the pure hopping model 
on a finite L x L periodic comb [6], with A^^. = L^. For Eq we have: 



Co{i,j) = Co{d)=^^*{i)^^{j)-^^L-'noiV2L)-' for L 

I — z 



• DO 



(28) 



Co(i, j) does not depend on the distance between the sites since the wave function of the ground state is constant 
along the backbone. In (28) L'^uq represents the filling of the state and {y/2L)~^ is the normalizing factor of the wave 
function. In the low energy hidden spectrum we have: 



(L-l)/4 

fe=i 



1 



L22cos(27rfc(i) 
{V2L) y/2l3tTT^k^ + 6 



E 



for L^oo (29) 



It should be noted that significant contribution to (29) arise only from states with energy close to Eq , similar to what 
happens in equation (18). Here the wave functions have macroscopical oscillations along the backbone and the filling 
of each state i^(\/2/3t7r^A;^ + 5)~^ is multiplied by an oscillating factor cos(27rA;rf). (\/2L)~^ is again the normalizing 
factor of the wave functions. In (Jq and (T+ we have (7^,-, (i,j) = Ccr^{i,j) = 0, since in this case tpj^ {i)'4'kU) ^s an 
oscillating factor with diverging frequency giving rise to decoherence effects. The large scale correlation function for 
the pure hopping model is then: 



Cii,j) = C{d) 



{s/2L) ' ^ V20tw^k'^ +Sj {y/2N;) ' ^ y/20tTT^k^ + d 



E 



2 cos(27rfcrf) 



Nr 



no 



E 



2 cos(27rfcd) 



for 



L 



(30) 



The two points correlation function diverges for L ^ cx), a well known property of localized condensates. Notice that 
each state of the condensate provides a difi^erent oscillating contribution to (30). In figure 5 we plot the correlation 
function C{i,j) as a function of the distance d. 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



FIG. 5. The correlation function given by (30) as a function of the macroscopic distance d. Above Tc, C{d) is zero, while at 
low temperature it depends on d. 



8 



We can also evaluate the correlation function between two sites of the same finger. Since all the states of macroscopic 
filling present the same behaviour along this direction (an exponential decay depending on the distance from the 
backbone), in this case there are no interference effects and one obtains an expression analogous to (27) and \j\ 
represent the distance of the sites from the backbone): 

(^/2L)p' V2Ptn^k^+S )~ iVm) tor r - L - oo (31) 

C{i,j) presents the diverging behaviour for L ^ oo typical of localized condensate. 



VI. THE GENERAL RESULT 



As we have seen, in an inhomogeneous structure the occurrence of BEC does not imply the macroscopic occupation 
of a single quantum state. In this last section we prove the general condition for a bosonic hopping model on a graph 
to present condensation only on the lowest energy state. In particular we will show that in a model presenting BEC 
at finite temperature (i.e with a low energy hidden spectra or with a spectral dimension greater than 2) one has 
condensation on the single state Eq if the spectrum p{E) satisfies the following spectral condition: 

where 

p"^{E) = ^Y!^^^-^'k) (33) 

^ k 

The symbol means that we are excluding from the sum the contribution given by the ground state. Equation (32) 
has a clear mathematical interpretation. Indeed as (4) is the necessary condition for evaluating the thermodynamic 
limit of the average value of any bounded function f{E}S) as / f{E)p{E)dE, (32) represents the condition to evaluate 
by an integral even unbounded functions diverging in Eq as {E — Eq)~^ . The general result can be obtained writing 
equation (6) as: 

/ = ^+/ '""IfL^'^'' -^^./ ^^.'E (34) 



- 1 ]e~„. =-'e"5 - 1 i-'ePB - 1 

where we set Eq = Q and e' is a generic energy between Eq and El the energy of the first excited state. Let us now 
take in (34) the thermodynamic limit r — > oo. For the first term limr^(x)(-2~^ — 1)~^ = no. For the second term we 
have: 

J ef^ - 1 - r^oo J \E\ ^ ' 

where we used the condition (32). Finally in the last term we let e' 0, obtaining: 

f p{E)dE 

no = f- J _ ^ = riTD (36) 

and this implies that only the lowest energy state has a macroscopic filling. The pure hopping model on the comb 
graph obviously does not satisfy (32). Notice that (32) is not a condition on p{E), the density of states of the infinite 
struc;tiires, but a prescription on how p^{E) tends to p(E). The single state occupation then is determined by the 
behaviour of the model on each large but finite mesoscopic scale. 



VII. CONCLUDING REMARKS 



In this paper we address the issue of Bose- Einstein condensation for a bosonic hopping model defined on the comb 
graph, a simple inhomogeneous discrete structure. Such model is expected to describe the low coupling limit for a 
comb-shaped network of Josephson junctions (or, possibly, of optical traps confining non-interacting atoms), a realistic 
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system which is currently being the subject of experimental research [4] . In particular we show that the macroscopic 
occupation of the ground state, which is one of the distinctive features of BEC on (homogeneous) Euclidean lattices, 
is not sufficient for a complete description of the phenomenon on the comb lattice. Indeed the inhomogeneity of the 
network forces the macroscopic occupation of a spectral region just above the ground state. This feature is strictly 
related to the presence of the so-called hidden regions in the spectrum of inhomogeneous structures such as the comb 
graph, and, as we discuss in Section V, it is the cause of an anomalous behaviour of the large distance two point 
correlation functions. These results are obtained by exactly evaluating the states filling on mesoscopic combs, where 
the energy spectrum is nearly continuous, but the levels are still distinguishable, being related to wave-functions which 
differ on a large scale. The conditions restoring the standard feature of the BEC, namely the macroscopic occupation 
of a single quantum state, are discussed as well. 
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